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We report a measurement of the B° — B° mixing parameter A md based on a 29.1 fb _1 sample of 
T(45') resonance decays collected by the Belle detector at the KEKB asymmetric e + e“ collider. We 
use events with a partially reconstructed B°(B ) —+ .D* =F 7r ± candidate and where the flavor of the 
accompanying B meson is identified by the charge of the lepton from a B°(B°) —> T^v decay. The 
proper-time difference between the two B mesons is determined from the distance between the two 
decay vertices. From a simultaneous fit to the proper-time distributions for the same-flavor (B°(B ), 
£±) and opposite-flavor (B°(B ), £ =F ) event samples, we measure the mass difference between the 
two mass eigenstates of the neutral B meson to be Amj= (0.509 ± 0.017 (stat)±0.020 (syst)) ps -1 . 
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I. INTRODUCTION 


After production, B° and B mesons evolve in time 
and mix into each other via the second-order weak inter¬ 
action box diagrams shown in Fig. |]. The mixing param- 
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FIG. 1: Standard Model “box diagrams” for the second-order 
weak B° — B mixing process. 


eter A rrid, which is the mass difference of the two neutral 
B mass eigenstates, is determined from the Feynman di¬ 
agrams shown in Fig. 0 ® to be 

Am d = ^f%mBm^rHS\V^V td \ 2 B B , ( 1 ) 

where m*, tub and mw are the t-quark, B° and W 
masses; Gf is the Fermi constant; rjt is a QCD correc¬ 
tion [||; S' is a function of to 2 jm\ v [||; fs is the de¬ 
cay constant of the B meson; and Bb is the B meson 
bag parameter |Q]. In principle, a measurement of the 
mixing parameter A rrid can be used to determine the 
magnitude of the CKM matrix element V t( i.\ . However, 
there are large theoretical uncertainties associated with 
the model dependence of fs and Bb- We report here 
a measurement of Amj that uses B mesons from the 
decays of T(4S) states produced by the KERB collider 
and recorded in the Belle detector. We determine the 
flavor of one B meson by partially reconstructing the de¬ 
cays B°{W) -> D* T 7r ^~; the flavor of the accompanying 
B meson is identified by the charge of the lepton from 
B°(B ) —> X^l±v decays. In the future, this technique 
can be extended to determine the linear combination of 
CKM angles 2</>i + p 3 ||. 

The T(45) decays to a B° — ~B° pair that is nearly at 
rest in the T(45 l ) center of mass system (CM). Since B 
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mesons are spin 0 mesons, angular momentum conserva¬ 
tion requires the two B mesons to be in an antisymmetric 
quantum state 

|a> = -±= (|H°(1))|#(2)) - |S°(1))|5 0 (2)>) , (2) 

where 1 and 2 indicate the opposite sides in the T(4S') 
decay plane. The form of Eq. (pi) guarantees that at 
any time the amplitude for either B°B° or B B states 
vanishes. Thus, if we determine the flavor and decay time 
ttag for one of the B mesons to decay into a final state 
ftag, then we can measure the time evolution of the other 
B at any time f as a function of the time difference At = 
t—ttag- For the two-state neutral B system, the following 
probability expressions hold to a good approximation [Q, 

|At| 

P OF (At) = jy—e t bo [1 + cos(Am d At)] , 

( 3 ) 

P SF (Af) = —F_e t bo [1 - cos(Am d Af)] , 

where the superscripts SF and OF denote events where 
the lepton-tagged and partially reconstructed B mesons 
have the same and opposite flavors, respectively. 

Since the B° and B are nearly at rest in the CM 
frame, At can be determined from the displacement be¬ 
tween the lepton-tagged and partially reconstructed B 
decay vertices 


At ~ ( z iec - z tag )/P^c = Az/P^fc, 


( 4 ) 


where the z axis is defined to be anti-parallel to the 
positron beam direction and the constant /3y = 0.425 
is the Lorentz boost of the e + e _ center of mass system 
at the KEKB collider. 

This paper is organized as follows: In Section [n], we 
describe the KEKB collider and the Belle detector. The 
partial reconstruction of B°(B ) —•> D* T ^ decays, the 
determination of the 6-flavor of the accompanying B me¬ 
son and the measurement of At are described in Sec¬ 


tion III. The likelihood fit to the measured At distribu¬ 
tions is described in Section 0- We present the results 
of the fit and studi es o f sources of systematic errors in 
Sections |V|. |VT and VII. The conclusions are presented 
in Section [VIII. 


II. EXPERIMENTAL APPARATUS 

KEKB is an asymmetric e + e _ collider 3 km in cir¬ 
cumference, which consists of 8 GeV e~ and 3.5 GeV e + 
storage rings and an injection linear accelerator. It has a 
single interaction point (IP) where the e + and e~ beams 
collide with a crossing angle of 22 mrad. The collider has 
reached a peak luminosity above 8xl0 33 cm _2 s _1 . Due 
to the energy asymmetry, the T(4S') resonance and its 
daughter B mesons are produced with a Lorentz boost 
of P'y =0.425. On average, the B mesons decay approxi¬ 
mately 200 /im from the T(4S') production point. 
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The Belle detector JtJ is a general-purpose large solid 
angle magnetic spectrometer surrounding the interaction 
point. Precision tracking and vertex measurements are 
provided by a silicon vertex detector (SVD) Q and a 
central drift chamber (CDC) jbj in a 1.5 T magnetic field 
parallel to the 2 -axis. The SVD consists of three layers of 
double-sided silicon strip detectors (DSSD) arranged in a 
barrel and covers 86% of the solid angle. The three layers 
at radii of 3.0, 4.5 and 6.0 cm surround the beam-pipe, a 
double-wall beryllium cylinder of 2.3 cm radius and 1 mm 
thickness. The strip pitches of each DSSD are 84 pm for 
the measurement of the 2 coordinate and 25 /mi for the 
measurement of the r — (f> coordinate. The CDC is a 
small-cell cylindrical drift chamber with 50 layers of an¬ 
ode wires including 18 layers of stereo wires. A low -Z gas 
mixture (He(50%)+C2H 6 (50%)) is used to minimize mul¬ 
tiple Coulomb scattering to ensure good momentum reso¬ 
lution, especially for low momentum particles. The CDC 
provides three-dimensional trajectories of charged parti¬ 
cles in the polar angle region 17° < 6 < 150° in the lab¬ 
oratory frame. The impact parameter resolution for re¬ 
constructed tracks is measured as a function of the track 
momentum p (measured in GeV/c) to be a xy = [19 © 
50/( p/3 sin 3 / 2 9)] /mi and a z = [36 © 42/(p/3 sin 5 / 2 6)] pm. 
The momentum resolution of the combined tracking sys¬ 
tem is a Pt /pt = (0.30//3©0.19pt)%, where p t is the trans¬ 
verse momentum in GeV/ c. 

The identification of charged pions and kaons uses 
three detector systems: the CDC measuring dE/dx , a 
set of time-of-flight counters (TOF) Q and a set of 
aerogel Cherenkov counters (ACC) [Tl 1 ]. The CDC mea¬ 
sures energy loss for charged particles with a resolution 
of a(dE/dx) = 6.9% for minimum-ionizing pions. The 
TOF consists of 128 plastic scintillators viewed on both 
ends by fine-mesh photo-multipliers that operate stably 
in the 1.5 T magnetic field. Their time resolution is 95 ps 
(rms), providing three standard deviation (3<r) 
separation below 1.0 GeV/c, and 2 a separation up to 
1.5 GeV/c. The ACC consists of 1188 aerogel blocks with 
refractive indices between 1.01 and 1.03 depending on 
the polar angle. Fine-mesh photo-multipliers detect the 
Cherenkov light. The effective number of photoelectrons 
is approximately 6 for /3 = 1 particles. Using the informa¬ 
tion from these three particle identification systems, the 
K /7r likelihood ratio P(K/n) = C(K)/(C(K) + C(n)) is 
calculated, where C(K) and C(i r) are kaon and pion like¬ 
lihoods 0. A selection with P(K/ir) > 0.6 retains about 
90% of the charged kaons with a charged pion misidenti- 
fication rate of about 6%. 

Photons and other neutral particles are reconstructed 
in a CsI(Tl) crystal calorimeter (ECL) [jl2| consisting of 
8736 crystal blocks, 16.2 radiation lengths (Xq) thick. 
Electron identification is based on a combination of 
dE/dx measurements in the CDC, the response of the 
ACC, the position and the shape of the electromagnetic 
shower, as well as the ratio of the cluster energy to 
the particle momentum [[b|. For the electron identi¬ 
fication requirement used in this analysis, the electron 


identification efficiency is determined from two-photon 
e + e _ —> e + e _ e + e _ processes to be more than 90% for 
Piab > 10 GeV/c. The hadron misidentification prob¬ 
ability, determined using tagged pions from inclusive 
K°s —> 7r + 7r decays, is below 0.5%. 

All the detectors mentioned above are inside a super¬ 
conducting solenoid of 1.7 m radius. The outermost 
spectrometer subsystem is a K/ and muon detector 
(KLM) |lj], that consists of 14 layers of iron (4.7 cm 
thick) absorber alternating with resistive plate counters 
(RPC). The KLM system covers polar angles between 
20 and 155 degrees. The efficiency of the muon identi¬ 
fication requirement used here, determined by using the 
two-photon process e + e _ —* e + e~ p + p~ and simulated 
muons embedded in BB candidate events, is greater than 
90% for tracks with pi a b > 1 GeV/c. The correspond¬ 
ing pion misidentification probability, determined using 
K°s —> 7r + 7r decays, is less than 2%. 

In our analysis, Monte Carlo (MC) events are gen¬ 
erated using the QQ event generator |l5|] and the re¬ 
sponse of the Belle detector is precisely simulated by a 
GEANT3-based program jlGj. The simulated events are 
then reconstructed and analyzed with the same proce¬ 
dure as is used for the real data. 


III. EVENT RECONSTRUCTION 

A. Data Sample 

We analyze a 29.1 fb -1 data sample recorded on the 
T(45) resonance. The data was taken from June 1999 to 
July 2001 and corresponds to about 3.13 x 10' BB pairs. 


B. BB Event Pre-selection 

In Belle, neutral B mesons can only be created via 
the process T(4S) —*• B°B . To suppress the non -bb 
background processes from QED, beam-gas and e + e _ —► 
t + t”, we select hadronic events using event multiplicity 
and total energy variables |nj . 

C. B —> D* 7r Decay Partial Reconstruction 

We now describe a partial reconstruction method for 
the decay chain B — > D*irf, D* —> Dir s , where 717 and 
7 r s designate a fast n and slow 7r, respectively. 1 This is a 
variation of a technique that was first developed by the 
CLEO collaboration |l^]. For the Belle experiment, we 
modify and apply this method to make a precise time 


1 Throughout the paper, the charge conjugate process is implied, 

e.g. B° —> D* _ 7r+ denotes also B° —> D* + tt~ etc. 
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dependent measurement of A nid- Unlike analyses that 
use fully reconstructed B —> . 0 * 77 / decays, we do not use 
any properties of the decay products of the D° meson 
in the decay D* —> Dtt s . According to a Monte Carlo 
simulation, this method yields an order of a magnitude 
more events than a full reconstruction method. 


The asterisks denote variables calculated in the D* rest 
frame. Calculated in this way, cos 9* can take values 
outside the physical region due to finite resolutions, or 
for background events. 

2. B°(B°) —► U* T tt* Event Selection 


1. Kinematics 


We consider five particles in the decay chain: 
B,D*,nf,D and n s . When reconstructed in this way, 
the system has 20 degrees of freedom. 

For partial reconstruction, only the 77/ and 7r s are used. 
The D candidate is not reconstructed. We can obtain 
constraints from 4-momentum conservation for both the 
decays B —> D*ttj and D* —> Dn s (8 constraints). The 
B 1 D *, D , 7 Tf and tt s masses from the PDG2000 compi¬ 
lation ||l9f provide further 5 constraints. In addition, we 
use the constraint that the B energy is the CM beam en¬ 
ergy of KEKB at the T(4S) divided by 2 (1 constraint). 
The measurements of the 7 7/ and 77 s 3-vectors provide 6 
constraints. In total there are 20 constraints, equal to 
the number of degrees of freedom of the system. Follow¬ 
ing previous analyses ([0| |2(J, we use two variables to 
measure the B —> 0*77 signal. The first is D° missing 
mass, 

Mz) miss = m 2 B + m^ f + m 2 a — 2E B E Vf — 2 E B E na 

T ‘2‘B 7T fE 7Ts T 2|p^j | \pttj | cos Obttj (5) 

+ 2\pb\ \Ptt s I cos - 21 p nf 11 p Wa | cos d 7If7Vs , 


where E, p, m are the energy, momentum and nominal 
masses for the B, D*, D , 717 and 7 r s mesons in the CM 
frame; 6 B n fM is the angle between the directions of mo¬ 
tion of the B and 7 r/( s ); and 0 nf7lTs is the angle between 
the 77 / and 7r s . In the CM frame, cos0b,t s ~ — cos d B Tr f ', 
the analysis |2(J shows that this approximation, and the 
relation 


COS 0b-k f 


-m 2 B - m 2 f + m 2 D , + 2 E B E Tf 
2\PB\\Pn f \ 


( 6 ) 


can be used to evaluate M£, miss without reconstructing 
the B meson’s flight direction. The second variable is the 
angle 9* a between the slow pion in the D* rest frame and 
the direction of motion of the D* in the CM frame. Using 
partial reconstruction, it is calculated using the relation, 


cos 61 = 


(3 d ,(E* d -E 1 s ) \p D \ 2 -\p„ 


m 


2jh,(3 D *mD* | Pi 


(7) 


where 7 £>• = = 

(e b - ^\pK f \ 2 /m D *, (3 D * = Vl - (1/7 d*) 

and \p D \ = \JE 2 d - m 2 D = yJ{E B - E„ f - E n J 2 - m 2 D ; 

E* d and El are the energy of D and 7 r s in the D* rest 
frame; p* a is the momentum of 717 in the D* rest frame. 


The B°(B°) D** 77"^" candidates are reconstructed 

with the following requirements. In the CM frame, we 
select nf candidates with momentum p nf in the range 
2.05 GeV/c < \p VJ \ < 2.45 GeV/c and an oppositely 
charged 7r s with momentum | p ns \ less than 0.45 GeV/c. 
We require dr < 0.05 cm, dz < 2 cm for the 77/ can¬ 
didate and dr < 0.2 cm, dz < 2 cm for the 7r s candi¬ 
date to suppress backgrounds from beam particles that 
interact with the residual gas of the vacuum system or 
spent-beam particles that strike the vacuum chamber 
wall. The variables dr and dz are the distances of clos¬ 
est approach of the track to the interaction point in the 
r — 4> and z planes, respectively. For better vertex res¬ 
olution, the SVD is required to provide at least 2 spa¬ 
tial points for a 77 / candidate track. The 77 / candidates 
are required to have electron and muon likelihood ra¬ 
tios less than 0.8 to suppress background from semilep- 
tonic B decays. Angular momentum conservation in the 
pseudoscalar to vector-pseudoscalar decay B° —> U*”7r + 
leads to a distribution proportional to cos 2 for the 
angle d* a . To enhance the signal to background ra¬ 
tio, we require 0.3 < |cos#*J < 1.05. We then select 
B° —> D*~tt + candidates with D° missing mass M B > miss 
greater than 1.85 GeV/c 2 and 0.3 < |cos0*J < 1.05. 
These two requirements define the signal region. 
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FIG. 2: cosd- Wf7ra — cosS-TrjU* cos£U s d* distributions in a 
simulation of fully reconstructed Monte Carlo events. 


Fig. shows that the signal peaks at zero in the vari- 
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able | cosfl-n^Tr, — cosO- nf D» cos9 Va D*\ For about 
4% of the events, more than one pair of opposite sign 
particles satisfy all the selection requirements. In these 
cases, we select the combination with the smallest value 
of | COS 0 — TT j^TT s COS 0 — j D * COS ^7T S Z)*|- @ — TTfK s 1 ^ — 7 T fD* 
and Ok s d* are the angles between the directions of —p-K fl 
Pd* and p JTa in the B —> D*Wf, D* —> Dn s decay, respec¬ 
tively. This quantity is zero when the direction of the 
7r s is a good approximation to the D* direction and is 
small when the plane defined by the pd * and p na nearly 
coincides with the plane defined by the pd* and —p nf • 


D. Flavor Tagging for neutral B meson 

The flavor of the signal B decay is obtained from the 
charge of the fast pion. The flavor of the accompany¬ 
ing B meson decay is determined from the charge of 
the primary lepton (e or p) from the semileptonic de¬ 
cay b —> cl~v. In addition to tagging the flavor of the B 
meson decay, there are two other major reasons for using 
a high-momentum lepton: the point of closest approach 
of the lepton track and the beam-line provides the loca¬ 
tion of the tagging side vertex; the requirement of a high 
momentum lepton in the event dramatically reduces the 
continuum background. 

We require the lepton to have a momentum in the CM 
frame of at least 1.1 GeV/c. We select leptons from well 
measured tracks by requiring dr < 0.05 cm, dz < 2 cm. 
We demand that the lepton have both r — (f> and z hits in 
the SVD. To reject secondary leptons from the decay of 
the unreconstructed D° mesons, we require the cosine of 
the angle between the lepton and ttj to be greater than 
—0.8. We also reject lepton tracks that, when combined 
with any other oppositely charged track in the event, 
have an invariant mass that is within 50 MeV of the J/^ 
mass. If more than one lepton in an event satisfies all 
of the above criteria, the highest momentum lepton is 
selected. 

However, not all lepton candidates are primary lep¬ 
tons. According to a MC study, the background is dom¬ 
inated by three sources. The first is secondary leptons 
from charm decays (“cascade leptons”) which come from 
the b —> c —> s£ + v decay chain and not directly from 
b —> ci~v like the “primary leptons”. The charges of 
the cascade and primary leptons are opposite and, thus, 
cascade leptons can bias the tagged flavor of B mesons 
and must be well understood. Two of the requirements 
described above discriminate against secondary leptons: 
the momentum requirement and the angular cut with re¬ 
spect to the fast n direction. The second category of in¬ 
correct tags is due to leptons from J/ij) and ip(2S) decays. 
Equal numbers of positively and negatively charged lep¬ 
tons are produced from J/ij) and ij){2S) decays. Hence, 
the charge of the observed lepton is uncorrelated with 
the B flavor. The third category is composed of hadronic 
tracks misidentified as leptons. Particle identification is 
applied to suppress this background. 


E. Vertex Reconstruction 

An analysis that relies on At information requires a 
measurement of At in Eq. (|j). We use the 2 difference 
between B decay vertices in the laboratory system to 
approximate the proper-time difference At by 

As = z nf - zi= t(A z C m + c(3At) ~ cfijAt , (8) 

where z nf and zi are the z positions of the fast pion and 
lepton tracks respectively, and c is the speed of light. 
The constant /3q = 0.425 is the Lorentz boost factor for 
the e + e _ center of mass system in the Belle experiment. 
A zcm is the z difference between B decay vertices in the 
CM frame. Here, we assume that the B mesons are at 
rest in the CM frame and thus A zqm — 0. Therefore we 
have 



( 9 ) 


The 2 positions are determined from the intersection 
of the 7Tf (or lepton track) with the profile of B decay 
vertices, which is estimated run-by-run from the pro¬ 
file of the interaction point (cr^ p ~ 110 pm, <jy P ~ 
5 pm. al p ~ 3500 pm) convolved with the average B 
flight length (30 pm in the CM frame). 


F. Signal 


Monte Carlo simulation shows that B°(B°) —> D*^p ± 
candidates are peaked around 1.865 GeV/c 2 in the 
D° missing mass distribution, the same location as 
B°(B°) —> D* T 7r ± decays. B°(B°) —> D*^p± decays can 
also be treated as signal for the measurement of B° — B 
mixing. Thus, the signal consists of B°(B ) —> D * =f 7t ± 
and B°(B ) —> D* zf p ± decays tagged by primary lep¬ 
tons. Henceforth, for simplicity we use the notation 
B°{B°) -> D**h ± to represent both B°(B°) D*^ tt* 

and B°(B°) —> D*^p ± decays. 

G. Backgrounds 

Some other decay modes, such as B° —> D**~n + and 
B° —> D*~l + vi, also peak in D° missing mass, and 
hence can fake the signal. We therefore divide the back¬ 
grounds into unpeaked and peaked categories. Unpeaked 
background is dominated by random combinations of nf 
and 7r s with primary leptons from B° and B± decays, 
and combinatorial background from continuum. We ver¬ 
ify in MC that the As shape of the unpeaked back¬ 
ground can be modeled by the D° missing mass side¬ 
band. The As shapes of unpeaked background thus are 
taken from the D° missing mass sideband. The peaked 
background is dominated by the following sources: signal 
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FIG. 3: D° missing mass (GeV) distribution for lepton tagged 
B° —^ D*~h + candidates. 


B decays ( B°(B °) 


D**w* and B°(B°) 


D* T p ± ) 

with secondary-lepton tags or fake-lepton tags; B° —► 
D**~ 7r + , B + —> D 7r + and B° —> I?* _ 7r + 7r 0 decays 
with primary-lepton tags, secondary-lepton tags or fake- 
lepton tags. The Az shapes of peaked background are 
determined from Monte Carlo simulation. The details of 
the background parameterization are described in Sec¬ 
tion m 


H. Signal Yields 


To obtain the signal yields, we apply a binned maxi¬ 
mum likelihood fit to the D° meson missing mass distri¬ 
bution. In the fit, the signal shape is determined from 
the B —> D*TTf Monte Carlo simulation. The shape of 
the unpeaked background is taken from wrong-sign com¬ 
binations, where the sign of the charges of the 7 r/ and 
7rs are the same. The MC simulation shows that the 
wrong sign combinations have a distribution consistent 
with the right sign unpeaked background. The peaked 
background shape comes from the BB Monte Carlo simu¬ 
lation. In the fit, the shape of the signal, the peaked back¬ 
ground and the unpeaked background are fixed. We also 
fix the peaked background normalization from the Monte 
Carlo simulation. We float the normalizations of the sig¬ 
nal and unpeaked background. By fitting the D° missing 
mass distributions shown in Figs |ii] and we obtain sig¬ 
nal yields of 3433 ± 81 (SF: 751 ± 41, OF: 2682 ± 70). 
The estimated backgrounds (unpeaked/peaked) in the 
signal region are 1466 ± 37 (SF: 243 ± 16/219 ± 15, OF: 
491 ± 18/513 ± 23) for M DmiBS > 1.85 GeV/c 2 . 

As a consistency check, we use another method to es¬ 
timate the signal yields. In Fig. ||, we fit the cos 9* s 
distribution with both signal and background shapes de¬ 
termined from MC. Here, we float both the signal and 


o 

> 

a> 

O 

fj 

o 

o 


in 

0) 


o 

> 

o 

O 


<N 

O 

o 

o 

in 

o 


c 

LU 



1.78 1.80 1.82 1.84 1.86 

D° missing mass (GeV/c 2 ) 


FIG. 4: D° missing mass (GeV/c 2 ) distribution (a) for oppo¬ 
site flavor final states and (b) for the same flavor final states. 


background normalizations. For the signal shape, we 
treat B° —> D*~ir + and B° —> D*~ p + events sepa¬ 
rately. The polarization of B° —> D*~p + is fixed from the 
CLEO measurement jin]]. We obtain B° —> H*~7r + and 
B° —> D*~p + yields of 2454 ± 90 and 1353 ± 96 events. 
Note that the yields include incorrectly tagged signal. 
According to a MC study, the probability of incorrect 
tagging is 9.5%. Thus, we conclude that the two meth¬ 
ods yield consistent results. Since the error is smaller 
for the yields obtained by fitting the D° missing mass 
distribution, we use the results from that method in the 
analysis. 


IV. MAXIMUM LIKELIHOOD FIT 


We extract the mixing frequency Ato^ by simultane¬ 
ously fitting the time evolution distribution of the SF and 










































FIG. 5: cos 9^ s distribution. The histogram from the fit and 
the contributions of signal and background are overlaid. 


OF B° —> D*~h + samples. The unbinned maximum like¬ 
lihood fitting method is applied to expressions containing 
Amd as a free parameter, which take into account both 
signal and background. 

A. PDF and Likelihood Function 

Here, we summarize the forms of signal and back¬ 
grounds used in the fitting. The likelihood is also es¬ 
tablished. 

For the signal, the probability density functions (PDF) 
of OF and SF events are given by 

/ OO 

P OF (Af')i? sig (A^ At')dAt' ,(10) 

-OO 

/ OO 

P SF (At')R sig (At — At')dAt' ,(11) 

-OO 

where f? s i g (Af) is the signal resolution function, which 
is parameterized by a triple-Gaussian distribution as fol¬ 
lows. 

I? si g(Af) = f 1 G(At;n 1 ,ai) + f 2 G(At-,fi 2 ,^2) 

+ (1 — h ~ f 2 )G(At; H 3 , cr 3 ) . (12) 

Here, G(t; n, a) is the Gaussian distribution, /i and f 2 
denote the fractions of the first and second Gaussian, 
respectively. For the background, the PDFs include sep¬ 
arate contributions for the parts that are peaked and 
unpeaked in the D° missing mass distribution. The time 
dependent parameterization includes prompt (zero life¬ 
time) and finite lifetime components as well as B° — B 


mixing. The time distribution of the unpeaked back¬ 
ground is determined from the D° missing mass sideband 
data. The functional forms for the background PDFs are 
given in the appendix. 

Using the PDFs described in Eqs. (lOi) ( |lT|) (|Al|) and 
(|A2|), the likelihood function can be written as 

£ = II ((! - /bkg)^? g F (Ati) + /bkg/bkgFbkg(Ati)) X 

i 

n ((! - /bk g )^g F (A tj ) + / bkg (l - / b ° k F g )F^ g (At,0)3) 
j 

where /bkg = ^bkg/(.Nbkg + N s i g ) is the background frac¬ 
tion. /bkg is the fraction of OF events in the background 
and calculated from TV^g/lVbkg- The normalizations are 
determined from the fit to the D° missing mass distribu¬ 
tion. 

B. Resolution Function 

We need to determine the detector resolution function 
to smear the theoretical probability expressions in Eq. [|. 
The signal resolution function R s i g (At) is the distribu¬ 
tion of the difference between the generated and recon¬ 
structed B° decay proper times. A good approximation 
to the resolution function can be obtained from the recon¬ 
structed As = Zi+ — Zi- distribution of J/ifi —> l + l~ de¬ 
cays in data, which were parameterized using the triple- 
Gaussian distribution in Eq. ©• In Monte Carlo simu¬ 
lation, Fig. H shows good agreement between the signal 
resolution function from a D*~n + + l^ ag ” events and the 
background subtracted As distribution from J/ip — > l + l~ 
decays. This indicates that using inclusive J/ip decays to 
model the tagged D*ir resolution function is justified. 



FIG. 6: Resolution functions for J/ip —> l + l decays and 
+ Itag” in MC. 
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J/ip —> l + l~ candidates were selected using similar re¬ 
quirements to those in the “B° —> D*~tt + + l// ag n selec¬ 
tion. The particle identification criteria are changed to 
select lepton pairs. Furthermore, no requirement is made 
on the angle between the tracks, since J/ip decays pro¬ 
duce lepton pairs with large opening angles . 2 In Fig. [i], 
the J/ip —> p+p _ mass distribution is fitted with the 
sum of a Gaussian and a second order polynomial. For 
the J/ip —> e + e^ case, the sum of a “Crystal Ball” func¬ 
tion | 2 ^| and a second order polynomial is used. 


we are able to extract the number of signal and back¬ 
ground events for both the e + e _ and p + p _ cases. To 
obtain all the parameters of the signal resolution func¬ 
tion, we apply a maximum likelihood fit to the overall 
Az distribution. In the Az fit, the background shape is 
determined from the upper sideband of dilepton mass. 
The normalization is obtained from a second order poly¬ 
nomial fit to the dilepton mass distribution. The back¬ 
ground shapes are also parameterized by triple-Gaussian 
functions. The parameters of the signal resolution func¬ 
tion (Eq. (p^)) are a± = (94±6) pm, a 2 = (227±18) /im, 
cr 3 = (736±98) pm, f x = 0.56±0.07and f 2 = 0.38±0.04. 
Here, the fit gives mean values, m =0, i = 1,2,3, for 
each Gaussian that are consistent with zero. 



FIG. 7: Invariant mass distributions for (a) J/ip —> p + p and 
(b) J/ip —> e + e _ . Fits are described in the text. 


By fitting the J/tp invariant mass distribution in Fig.|f], 



■ 1600 . - 800 . 0 . 800 . 1600 . 

Az(|im) 


FIG. 8 : Signal resolution function fit with a triple-Gaussian 
for J/ip —> l + l~ decays in the data. Plot (a) shows the fit on 
a linear scale and plot (b) shows the fit on a semilogarithmic 
scale. 


Similarly, after applying the same technique to MC 
samples of J/ip —> l + l~ decays, we find the parameters 
of the MC resolution function i?Mc(At) to be 01 = (79 ± 
4) pm, 02 = (203 ± 13) pm, 03 = (625 ± 74) pm, fi = 


2 


Removing this requirement is found not to affect the measured 
resolution function. 
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0.49 ± 0.08 and /2 = 0.43 ± 0.04. We also fix the mean 
values of /t* = 0, i = 1,2,3 to zero. 

From the J/ij} —> l + l~ sample (see Fig. |i)|), we observe 
that there is a discrepancy between the distributions in 
data and MC simulation. To account for this difference, 
we convolve the resolution function in MC with a Gaus¬ 
sian of 50 /im width to reproduce the resolution in the 
data. This additional smearing is applied to the MC sim¬ 
ulations of peaked backgrounds. This procedure was ver¬ 
ified by comparing the Az distributions of the D° missing 
mass sideband in data and MC. 



-1500 -1000 -500 0 500 1000 1500 


Az (gm) 


FIG. 9: A z distribution for J/%j} —> l + l decays in data and 
Monte Carlo without (a) and with (b) the extra 50 /im smear¬ 
ing. 


For the peaked background case, we use the resolution 
function for J/ip —> l + l~ decays in MC with an additional 
50 /m i smearing (i?^ k (At) = i?Mc(At) ® G(At;/z = 
0,cr = 50 /im)). With the additional 50 fi m smear¬ 
ing, we extract the resolution function of the MC sim¬ 
ulation. We find that the mean value of the second 
Gaussian of the resolution function for the same fla¬ 


vor peaked background events is inconsistent with zero, 
H 2 = (81 ± 45) /im. The corresponding mean value for 
the opposite flavor peaked background is consistent with 
zero. Thus an offset is included only in the same flavor 
peaked background shape. 


C. Background Distributions 


The parame teri zation s of all background shapes are 
given in Eqs. (Al) and (|aJ). Figs. [To] and [ll| show the 
A z distributions for unpeaked and peaked backgrounds, 
respectively. The unbinned maximum likelihood fitting 
method is applied to extract the background parameters. 



FIG. 10: Semilogarithmic plot of the Az distributions from 
the D° missing mass sideband for unpeaked backgrounds in 
(a) opposite-flavor and (b) same-flavor events. 


We use the Az distribution from the D° missing mass 
sideband to reproduce the unpeaked background shape. 

For the Az shapes of the peaked background, we use 
the Monte Carlo simulation. We find that the probability 
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FIG. 11: Semilogarithmic plot of the A z distributions from 
MC simulation for peaked backgrounds for (a) opposite-flavor 
and (b) same-flavor events. 


• The B° lifetime, Tgo = (1.548 ± 0.032) ps, is fixed 
to the PDG2000 value JuJ . 

• The signal and unpeaked background resolution 
function are determined from J/ip —» l + l~ de¬ 
cays in data. The parameters of Eq. (0) are: 
(7 1 = (94 ± 6 ) fjm, <j-2 = (227 ± 18) /rm, 173 = 
(736±98) /tm, /1 = 0.56±0.07 and f 2 = 0.38±0.04. 
In the fit, all the parameters are fixed to their cen¬ 
tral values. 

• The peaked background resolution function is de¬ 
termined from J/ip —> l + l~ decays in MC con¬ 
volved with an additional 50 micron smearing term. 
The parameters of Eq. (^2|) are: <7i = (93 ± 4) fin 1 , 
(72 = (209 ± 13) Atm, a 3 = (627 ± 74) Atm, f x = 
0.49 ± 0.08 and / 2 = 0.43 ± 0.04. In the fit, 
all the parameters are fixed to their central val¬ 
ues. Note that the offset of the second Gaussian 
is At 2 = (81 ± 45) Atm for the same flavor peaked 
background. 

• The parameters of the unpeaked background 
shapes are fixed to values determined from fitting 
the D° missing mass sideband (see Fig. |lo|). The 
parameters for opposite flavor events are: fo = 
0.37±0.06, 7bk g = (1.98±0.17) ps. For same flavor 
events: fo = 0.25 ± 0.06, Tbkg = (1.41 ± 0.10) ps. 

• The parameters of the peaked background shape 
are fixed to values determined from fitting the 
Monte Carlo simulation in Fig. [□]. Here, A rrid is 
fixed to 0.509 ps -1 . The parameters for opposite 
flavor events are: fo = 0.07 ± 0.07, Tbkg = (1.34 ± 
0.12) ps. The fit gives f\ = 0. For same flavor 
events: fo = 0.10 ± 0.04, /1 = 0.40 ± 0.05, Tbkg = 
(1.84 ±0.26) ps. 


given by Eq. (|A5|) for the opposite flavor sample gives a 
good fit even without including the mixing term. The 
fraction of mixing, / 1 , is then set to zero and fixed for the 
fit to the opposite flavor sample. Therefore, the floating 
parameters are fo and Tbkg- For the same flavor events, 
inclusion of the mixing term is necessary to obtain a good 
fit and the floating parameters are fo, f± and Tbkg- The 
mixing frequency in the background is determined from 
the fit. 


V. FITTING RESULTS 

We now discuss the result of the fit. In the final fit, 
A md is the only free parameter. The parameters of the 
resolution functions and background shapes are fixed to 
their central values. The uncertainties in these parame¬ 
ters are included in the systematic error. The parameters 
used in the fit are the following: 


• The fractions / bkg = 0.30T0.01, f$£ = 0.69T0.06, 

/bkg~ Unpeak = °- 49 ± 0.03 and /bk g ~ Unpeak = °- 53 ± 
0.04 are determined by fitting the D° missing mass 
distribution in Fig. [|. 

Using the parameters determined above, Ato^ is ex¬ 
tracted from the fit. We find 


A md = (0.509 ± 0.017 (stat)) ps 


(14) 


In Fig. [T| we show the A 2 distributions for SF and OF 
data together with the curves from the fit. To display 
the charge asymmetry (see Fig. 13) between SF and OF 
events, we use 


A(Az) 


N OF (As) -N sf (A^) 
N of (A 2 ) ±N SF (Az) ’ 


(15) 


where N(Az) is the yield of the signal candidates as a 
function of A 2 . 
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Az (gm) 

FIG. 12: A z distributions for (a) opposite-flavor and (b) 
same-flavor events in data. The curve from the fit and the 
contributions of signal, unpeaked and peaked backgrounds 
are overlaid. 


VI. VALIDATION CHECK 



FIG. 13: Distribution of the asymmetry, A(Az), as a function 
of A z for data with the fit curve overlaid. 


TABLE I: Summary of signal MC bias test. 


A rrid (ps J ) for MC generation 

A rrid (ps x ) from fit yields 

0.442 

0.448 ±0.011 

0.472 

0.470 ±0.012 

0.502 

0.492 ±0.011 

0.532 

0.537 ±0.012 


We have also examined the variation of the fit results 
as the values of the cuts on the lepton momentum p* 
in the CM frame and the fiducial angle requirement for 
the lepton and nf. The results are shown in Figs 0(c), 
(d). We find that all the variations are consistent with 
statistical fluctuations. 

Furthermore, we performed fits separately for SF and 
OF events. Table [n| shows the fit results for SF events 
and OF events, respectively. 


We also perform a series of fits to signal Monte Carlo 
using the same procedure used in the analysis of the data. 
We generate four signal MC samples with different A 
values. Table | summarizes the fitting results and shows 
the consistency between the input and extracted A md- 
We conclude there is no significant fit bias. 

We also check for run dependence of the Amj extrac¬ 
tion. The value of A rrid is shown for four different exper¬ 
imental data taking periods (referred to as experiments 
7, 9, 11, 13). The results are shown in Fig. 0(a). 

To check the dependence of the A rrid fitting on the 
continuum background level, we also calculate Ato^ as 
a function of the cut on the normalized second Fox- 
Wolfram moment |23j i ?2 = H 2 /H 0 . The results are 
shown in Fig. |l4|(b). 


TABLE II: Summary of A nid fit results for SF (OF) events. 
Only statistical errors are shown. 


Sample 

Events 

A m d (ps -1 ) 

SF 

1213 

A rrid = 0.511 ±0.035 

OF 

3686 

A m d = 0.504 ± 0.019 


We also fit e tagged events and n tagged events sepa¬ 
rately (see Table |I|) . 

To check the sensitivity of the result to tails of the 
vertex resolution function, we vary the A z range of the 
fit. The results are shown in Table IV and are consistent 
with the primary result. 
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FIG. 14: A rrid as function of (a) experiment number, (b) R 2 cut value, (c) cut on the lepton momentum in the CM frame, (d) 
lepton fiducial angle requirement. 


TABLE III: Summary of Afit results for e and p tagged 
events. Only statistical errors are shown. 


Sample 

Events 

Arrid (ps -1 ) 

e tagged 

2324 

Arrid = 0.510 ± 0.025 

p tagged 

2575 

Arrid = 0.503 ± 0.024 


TABLE IV: Summary of Arrid fit results for different Az win¬ 
dow. Only statistical errors are shown. 


max(|Az|) (pm) 

Arrid (ps -1 ) 

1700 

0.509 ±0.017 

1275 

0.507 ±0.017 

875 

0.511 ±0.019 

425 

0.520 ± 0.024 
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VII. SYSTEMATIC UNCERTAINTIES 

Various possible sources of systematic uncertainties 
were investigated. In Table [v| we list the contributions 
to the systematic errors in Athat were estimated by 
varying the relevant parameters by one standard devia¬ 
tion. 


TABLE V: Contributions to the systematic error. 


Source 

Errors (ps x ) 

Background fraction 

±0.014 

Signal resolution function 

±0.012 

Background shape 

±0.005 

B° lifetime (1.548 ± 0.032 ps) 

±0.005 

Detector resolution (50 ± 18 /zm) 

±0.002 


Total ±0.020 


• Background fraction_ 

As described in Section IV A, the background parts 
of the likelihood function are weighted by /bkg, 

fZ~ unpeak and fZ~ unpeak in Ec i- ©> 

which are determined from a fit to the D° missing 


mass distribution in data. The corresponding sys¬ 
tematic error, ±0.014, is estimated by varying the 
parameters of these background fractions by ±ler. 


• Signal resolution function 

The fitted parameters of the Az distributions for 
J/iJ) —» l + l~ events are varied according to the lim¬ 
ited data statistics by ±1 <t. The difference from 
the standard fit is taken as a measure of system¬ 
atic error for each parameter. The corresponding 
uncertainty for each parameter is added in quadra¬ 
ture in order to estimate a systematic uncertainty 
in A md of ±0.012 ps -1 . 


• Background shape 

For both peaked and unpeaked backgrounds, the 
Az distributions are determined from a method de¬ 
scribed in Section IV C| . Varying all parameters 
by ±lcr, we take the differences, ±0.004, from the 
standard fit as the systematic error. 


• B° lifetime 

The value of the B° lifetime was varied accord¬ 
ing to the uncertainty of the PDG2000 jlij value. 
Changes in A/tt^ of ±0.005 are observed. 


• Detector resol ution 

In Section IV C , we mentioned that the MC back¬ 
ground distributions are smeared to correct for the 
difference in the vertex resolution between the MC 
prediction and the data. The smearing function 
was a Gaussian with er = 50 /irn. We varied the 
amount of smearing by its uncertainty (18 //in) J24J] 
and repeated the fit. We obtain a difference in A nid 
of ±0.002 


The total systematic error, obtained by summing all 
errors from the different sources in quadrature, is 
±0.020 ps" 1 . 


VIII. CONCLUSION 

Using 29.1 fb -1 of data collected with the Belle detec¬ 
tor at the T(4S I ), we have measured the B° — B mix¬ 
ing frequency Am-d in B°(B°) —> D * T 7r ± decay with a 
partial reconstruction technique. The data were accumu¬ 
lated between January 2000 and July 2001. 

The asymmetric e + e _ beam energies of the KERB col¬ 
lider allows for the extraction of the time evolution of the 
B meson wave function from precise measurements of the 
decay vertex positions. The data are separated into SF 
and OF samples. The B° decay vertex resolution for 
the signal was determined from the Az distribution of 
J/ij) —> l + l~ decays that occur in the same data sam¬ 
ple. The backgrounds are divided into two components, 
peaked and unpeaked backgrounds. The Ato^ value, ob¬ 
tained by simultaneously fitting the SF and OF time dis¬ 
tributions, is 

Amd = (0.509 ± 0.017(stat) ± 0.020(syst)) ps -1 . 

This is the first direct measurement of Amd using the 
technique of partial reconstruction of B°(B ) —> D * =F 7r ± 
decays. This measurement is statistically uncorrelated 
with all other reported experimental results for Amd- 
It is also almost systematically independent from all 
other measurements. The systematic uncertainty is dom¬ 
inated by the background fractions and the signal reso¬ 
lution function; all quantities are measured experimen¬ 
tally except for the peaked background. This measure¬ 
ment agrees with the world average value of Amd = 
(0.472 ±0.017) ps -1 |jl9), and serves as a validation of the 
technique that will be used in the future for the measure¬ 
ment of the CP violation parameter sin(2</>i + ^>3) [ po| . 
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APPENDIX A: BACKGROUND PDFS 


For the background, the PDFs of OF and SF events 
are parameterized by 


r c 

^b°kg(A/) = / 

J —< 


•OF—unpeakpOF—unpeak/ a ,/\ 
bkg r bkg ) 


x ic;r ak (Af - At') 


->unpeak / 
g ^ 

/•OF-unpeak^pOF—peak ( 
' bk 
-»peak / 


, /-i _ /-UF-unpeaKxpUF-peaK /a ,/\ 

+ V 1 _ J bkg / r bkg ) 


X i? p ° ak (At - At'))dAt' , (Al) 


^bkg(Al) = 


x A“ peak (Af - At') 

/-i _ /*SF—unpeak\pSF—peak /a ./\ 

V 1 J bkg / r bkg \^ z ) 

X RllfiAt - At'))dAt' . (A2) 


Here, f^ g ~ unpeak = and / b ° k F " unpeak ee 


’bkg 


bkg 


' bkg 


Ar OF-u npeak / iv O k F 

are the fractions of unpeaked back- 


’ bkg 


ground in the SF and OF events. N is the number of 
events, the subscripts bkg and sig denote the background 
and signal components. f? bk ak (At) and f? bk peak (Af) 
are resolution functions for peaked and unpeaked back¬ 
ground and are parameterized by triple-Gaussian distri¬ 
butions, which have the same form as the signal resolu¬ 
tion function R s i g (At) in Eq. (|lj). In the fit, we use the 
signal resolution function for f? bk pcak (Af). The peaked 
background resolution function i? bk ak (At), described in 
Section IVB, is determined from J/ip —> l + l~ decays in 
MC with additional smearing. 


In the unpeaked background PDF, P 
and P bkE / unpeak (At) are defined by 


OF—unpeak 
bkg 


(At) 


P° k F g - unpeak (At) = f° F ~ unpeak S(At) 

~h(l yOF—unpeak^ 


2t, 


bkg 


-e Tbk « (A3) 


Pbkg" Unpeak ( Ai ) = /o SF-unpeak <5(Af) 


_j_^ ySF—unpeak^ 


1 |At[ 

e Tbk « (A4) 


2Tbki 


In the peaked background PDF, P b] ^ peak (At) and 
P bkg ~ peak (Af) are defined by 


p ° F g - peak (At) = / 0 OF - peak <5(Af) 


| ^ yOF—peakj^ yOF—peak^ 


l 


■ /-| />OF—peak\ /»OF—peak 

— Jo )J i 


2rbk g 

1 + A m%T& 


|Atj 

- e ’'’bkg 


d 1 bkg 


4Ttkg + 2Am^r bkg 
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|At] 

x e Tbkg [1 + cos(AmdAt)] , (A5) 

pSF- p ea k ( A £) = /p F_peak (5(At) 

1 I At| 

+(1 - /o SF_peak )(l - yS F- P eak )—I— e _7 bkl 

2rbk g 


■ /-i />SF—peakx /»SF—peak 

“•"V 1 — Jo )J 1 


1 + Ayn^kg 

2Am d T 'bkg 


|At| 

x e Xbkg [1 — cos(Am<jAi)] . (A6) 


In the above equations, /o is the prompt lifetime fraction 
and (1—/o)/i is the mixing fraction. 7bk g is the lifetime of 
the background component that does not originate from 
a prompt source. A rrid contributes to the peaked back¬ 
ground since the peaked background is dominated by B° 
decays (2(J. Thus, mixing terms are included in the PDFs 
for the peaked background. In the unpeaked background 
expression, we set the mixing term to zero since a pa¬ 
rameterization without mixing reproduces the data well. 




